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Abstract. Only 235 entries were processed as astrometric bi- 
naries with orbits in the Hipparcos and Tycho Catalogue (ESA 
1997). However, the Intermediate Astrometric Data (IAD) and 
Transit Data (TD) made available by ESA make it possible to 
re-process the stars that turned out to be spectroscopic bina- 
ries after the completion of the Catalogue. This paper illus- 
trates how TD and IAD may be used in conjunction with the 
orbital parameters of spectroscopic binaries to derive astromet- 
ric parameters. The five astrometric and four orbital parameters 
(not already known from the spectroscopic orbit) are derived 
by minimizing an objective function (x 2 ) with an algorithm of 
global optimization. This code has been applied to 8 1 systems 
for which spectroscopic orbits became available recently and 
that belong to various families of chemically-peculiar red gi- 
ants (namely, dwarf barium stars, strong and mild barium stars, 
CH stars, and Tc-poor S stars). Among these 81 systems, 23 
yield reliable astrometric orbits. These 23 systems make it pos- 
sible to evaluate on real data the so-called 'cosmic error' de- 



scribed by Wielen et al. (1997), namely the fact that an un- 
recognized orbital motion introduces a systematic error on the 
proper motion. Comparison of the proper motion from the Hip- 
parcos catalogue with that re-derived in the present work indi- 
cates that the former are indeed far off the present value for 
binaries with periods in the range 3 to ~ 8 years. Hipparcos 
parallaxes of unrecognized spectroscopic binaries turn out to 
be reliable, except for systems with periods close to 1 year, as 
expected. Finally, we show that, even when a complete orbital 
revolution was observed by Hipparcos, the inclination is unfor- 
tunately seldom precise. 



Key words: Methods: data analysis - Astrometry - binaries: 
spectroscopic - Stars: late-type - Stars: fundamental parame- 
ters - Stars: barium, S, CH 



Send offprint requests to: pourbaix@astro.ulb.ac.be 

* Based on observations from the Hipparcos astrometric satellite 
operated by the European Space Agency (ESA 1997) 

** ESA-PRODEX contract 133 1 8/98/NL/VJ 
*** Scientific Collaborator, F.N.R.S., Belgium 

* Research Associate, F.N.R.S., Belgium 



1. Introduction 

The chemical anomalies observed in several classes of late- 
type stars are due to mass transfer across a binary system. This 
scenario, first suggested by McClure (1983)| , holds for dwarf 
and giant barium stars (Pop. I G and K stars with overabun- 
dance of carbon and heavy elements like Sr and Ba produced by 
the s-process of nucleosynthesis), CH stars (the Pop. II analogs 
of Ba stars), and extrinsic S stars (late-type giants character- 
ized by ZrO bands and no Tc lines, an element with no stable 
isotopes). The set of spectroscopic orbits available for those 
chemically-peculiar red stars (CPRS) has been considerably 
enlarged thanks to a decade-long radial- velocity monitoring. 
Such orbits are no w available for 63 gi ant barium stars and 19 
extrinsic S stars (lorissen et al. 1998), 21 dwarf barium and 



subgiant CH stars (McClure 1997;North , priv. comm.) and 8 
CH stars ([McClure & Woodsworth 1990b. Most of them were 



not available to the Hipparcos reduction consortia. The astro- 
metric parameters (parallaxes and proper motions) published 
in the Hipparcos catalogue (ESA 1997) were therefore derived 
from a single star solution (the so-called '5 -parameter model'). 
Since many of these CPRS have orbital periods in the range 1 
- 5 yr, their unaccounted orbital motion may seriously confuse 
the parallax and proper motion as determined by the Hipparcos 
consortia. Thus, unrecognized orbital motions might introduce 
a so-called 'cosmic error' (of a systematic nature) in the proper 



motion (Wielen 1997; Wielen et al. 1999). Parallaxes may also 



possibly be in error, especially for binaries with orbital periods 
close to 1 yr. 

The conclusions of previous studies inferring the kinemat- 



ics and absolute magnitudes of CPRS ( 


Bergeat&Knapikl997 


Mennessier et al. 1997 




Van Eck et al. 1998h from the data 



of the Hipparcos catalogue may possibly be affected by these 
systematic errors. It is therefore of importance to re-evaluate 
the astrometric parameters of CPRS with a model correctly ac- 
counting for the orbital motion (the so-called '12-parameter 
model', reducing to 9 parameters if the spectroscopic orbital 
parameters P, T and e are fixed), and to compare these values 
with those from the Hipparcos catalogue to identify the stars 
most affected by these systematic errors. 
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This comparison should allow to assess the reliability of 
the kinematic properties and absolute magnitudes of CPRS 
formerly derived from the Hipparcos-catalogue values. As a 
by-product, the astrometric orbit has been obtained in many 
cases, leading to an estimate of the masses through the ratio 
M 2 /[l + (Afi/M 2 )] 2 = /(M)/sin 3 i (derived by combin- 
ing the orbital inclination from the astrometric orbit with the 
mass function f(M) from the spectroscopic orbit when only 
one spectrum is observable). 

From the operational point of view, the numerical re- 
processing of the raw Hipparcos data (Intermediate Astromet- 
ric Data or Transit Data; see Sect. 2.1) presented in this paper 
is innovative on the following points: (i) a global optimization 
technique is used to minimize the objective function in the 12- 
parameter (or, for spectroscopic binaries, 9-parameter) space, 
and (ii) a change of independent variable allows to derive al- 
ways positive parallaxes, contrarily to the situation prevailing 
in the Hipparcos catalogue. The tool developed here for CPRS 
may in the future be applied to any other binary star to re- 
process the Hipparcos data when new orbits appear. The sec- 
ond paper of this series will be devoted to the re-processing 
of spectroscopic binaries involving a giant component (Boffin 



et al. 1993). Another is devoted to Procyon ( Girard et al. 2000) 



We encourage interested readers to communicate to the authors 
other spectroscopic binaries whose Hipparcos data would need 
to be reprocessed. 



2. Numerical methods 

In order to allow some re-processing of the Hipparcos data, the 
Hipparcos consortia provided the users with two kinds of data: 
the Transit Data (TD) and the Intermediate Astrometric Data 



(IAD). TD (Quist & Lindegren 1999) are a by-product of the 



analysis performed by the NDAC consortium, and merge as- 
trometric and photometric data. They basically provide the sig- 
nal as modulated by the grid in front of the detector. They are 
therefore especially well suited to the analysis of visual double 
or multiple systems, since the separation, position angle and 
magnitudes of the various components may in principle be ex- 
tracted from the TD. For that reason, TD are not provided for all 
Hipparcos entries, but only for those stars known or suspected 
of being double or multiple systems. IAD are lower level data 
which provide astrometric information only (the star's abscissa 
along a reference great circle, and t he pole of that circle on 
the sky; van Leeuwen & Evans 1998), but for every Hipparcos 
entry. In the following, we describe how IAD and TD (when 
available) have been used to derive the astrometric parameters 
of the CPRS that have a known spectroscopic orbit. When TD 
are available, the astrometric parameters may be derived in- 
dependently from the two data sets, and allow an interesting 
internal consistency check. In both cases, the method proceeds 
along the following steps: 

- definition of an objective function to minimize 

- computation of the objective function from the astrometric 
parameters 



- global minimization of the objective function with the sim- 
ulated annealing algorithm 

- local minimization of the objective function (with a quasi- 
Newton method) 

- comparison of the values of the objective function for 
the Hipparcos 5-parameter model and for our 9-parameter 
model, and application of an F-test to evaluate whether the 
9-parameter model is significantly better than the Hippar- 
cos solution 

- check of the astrophysical consistency of the 9-parameter 
solution from some a priori knowledge of the properties of 
the system (i.e., component masses) 

These steps are discussed in turn in the remaining of this sec- 
tion, first for IAD, and then for TD. 



2.1. Intermediate Astrometric Data 
2.1.1. Objective function 

The first step in any data-fitting problem is to set up an ob- 
jective function in order to compare different solutions (cor- 
responding to different values of the model parameters). This 
function is usually constructed in such a way that its lowest 
value corresponds to the best solution. The problem then re- 
duces to minimizing that function. The IAD already include 
some corrections (e.g. aberration and satellite-attitude correc- 
tions) and they are no longer stricto sensu raw data. They pro- 
vide the abscissa residuals along a reference great circle for 
the Hipparcos 5-parameter solution (i.e. ao, Sq, vd, /J, a , us, re- 
spectively: the position in right ascension and declination in 
the equinox 2000.0 at the epoch 1991.25, the parallax, the 
proper motions in right ascension and declination; following 
the practice of the Hipparcos catalogue, aj — ao cos 5 and 
/i* = [i a cos 5 will be used instead of ao and fi a ). By slightly 
mod ifying these 5 par ameters and by adding 7 orbital parame- 
ters ( Binnendijk 196C ), the aim is to further reduce the abscissa 
residuals Av below the values obtained from the Hipparcos 5- 
parameter model. In the case of the Hipparcos Intermediate As- 
trometric Data, the objective function is the x 2 expressed by 
Eq. 17.12 of Vol.3 in ^SA(1997)| 



M 



(Aw _^|!L Apfc )V-i (At; _^|_L Apfe)) (1) 



k=i 



k=l 



where the superscript t means transposed, Av are the abscissa 
residuals provided by the IAD file and corresponding to the 
Hipparcos 5-parameter solution, and is the partial deriva- 
tive of the abscissa with respect to the fc-th parameter express- 
ing how the abscissa residual varies when a correction Apk is 
applied to the value of the fc-th parameter with respect to the 
Hipparcos solution. M is the number of parameters retained in 
the solution, and V^ 1 is the inverse of the co variance matrix of 
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the observations given by: 



V 



(Vx 

v 2 
V o ••• 



'■. 

14/ 



with 



(2) 



If an observation j was processed by the two consortia (FAST 
and NDAC), the residuals obtained by the two reduction con- 
sortia are correlated and Vj is the 2x2 variance-covariance 
matrix for measurement j. On the other hand, when only one 
consortium processed observation j, Vj reduces to one num- 
ber^] (the estimated uncertainty aj of the measurement). Av, 
together with J^- (k = 1 to 5, with p\ = a^,P2 = S,p3 = 
va,p4 = fJ, a *,P5 = us), the original astrometric parameters 
as well as p, ap. and apf. are given in the IAD file. In order 
to evaluate \ 2 (Eq. [I]) for an orbital model, expressions for the 
partial derivatives of v with respect to the orbital parameters 
are required. They can be expressed as a function of the partial 
derivatives of v with respect to and So as follows (Hipparcos 
catalogue, Vol. 3, Eq. 17.15): 



dv 
do 



dv <9£ 
da% do 



dv dr] 
d5o do 



where o is any orbital parameter and 

£ = a* + p a , (t - t ) + RP a w + y, 
r] = So + ps{t — t ) + RP$vj + x. 



(3) 



In the above expression, £ and T) represent the Cartesian coor- 
dinates of the observed component on the plane tangent to the 
sky at the position (ceo, So). They combine the displacements 
due to the proper motion, the orbital motion and the parallax. 
R is the radius vector of the Earth's orbit in A.U. at time t. P a 



and Ps, the parallax factors, are given by Binnendijk (1960) 



P a = cos e cos a sin — sin a sin 
Ps = (sin e cos 8 — cos e sin a sin S) sin 
— cos a sin S cos 

where and e are respectively the longitude of the Sun and the 
obliquity of the ecliptic, both at time t.The variables x and y 
describe the apparent orbit (i.e., projected on the plane orthog- 
onal to the line of sight) of the observable component around 
the center of mass of the system. They are usually expressed 

1 In order to have a unique expression for \ 2 regardless of the 
dataset available (FAST, NDAC or both), we always write Vj as a 
2 by 2 matrix. When, for whatever reason, a measurement from only 
one consortium is used, that observation is duplicated and its weight 
divided by two (Vj = diag(2<jj , 2a j )). There is no change on the 
actual value of \ 2 but its expression is easier (though not faster!) to 
evaluate. 



in terms of the Thiele-Innes constants A, B, F, G of the photo- 
centric orbit as ( Heintz 1978| ) 



x = AX + FY 
y = BX + GY 



(4) 



with 

X = cos E — e 

Y = yjl-e 2 sinE, 

where E is the eccentric anomaly and (X, Y) are the coordi- 
nates in the true orbit. It should be noted that [x, y) and (X, Y) 
are referred to a Cartesian system with x pointing towards the 
North ( Heintz 1978| ), contrary to (£, rf) where £ points towards 
increasing right ascensions. For all the systems considered in 
the present paper, the magnitude difference between the two 
components is larger than the Hipparcos detection threshold 
(since the companion to CPRS is a cool white dwarf; Sect. || 
and Jorissen et al. 1998). Thus, one can safely assume that there 
is no light coming from the secondary. Hence, the orbit of the 
photocenter of the primary component is the same as the ab- 
solute orbit of the primary around the center of mass of the 
system. 

At this point, it is important to realize that the Hipparcos 5- 
parameter solution, that is used as a starting point for the new 
12-parameter solution, is in fact equivalent to a 12-parameter 
solution where the semi-major axis of the orbit (do) is null. It 
is therefore enough to consider in Eq. |l]the correction term rel- 
ative to the semi-major axis of the orbit (and in fact Aa = a , 
since the initial value of ao is null). The other orbital param- 
eters i, u, fl, e, P and T enter Eq. [I] only through the partial 
derivatives -J^- and -S^- [equal, respectively, to -J^ and -S^- 

according to Eq. |J entering These other orbital parame- 
ters do not require explicit correction terms in Eq. [I] (and their 
starting values would be ill-defined anyway). Hence, the orbital 
solution is the one minimizing 



where 



dv 



3 = Av - V j—Apk - ( 
^ dp k 



k=l 



dv <9£ ^ dv drj 
1 dpi dao dp2 dao 



(5) 



)ao (6) 



is a vector of dimension N (N is the number of observations), 
and £ and rj are functions of ao, i, ui, O, e, P and T, and J^- = 

r? 2 - = M^-. Thus, the 12 parameters enter the evaluation 
of x 2 although there are only six correction terms subtracted 
from Av. 

Our experience has shown that, except for very special 
cases (i.e., parallaxes and semi-major axes larger than about 
20 mas, orbital periods significantly different from 1 year but 
smaller than 3 years; one example is HIP 50805 in Table [l]), as- 
trometric orbits could not be derived from the IAD without an a 
priori knowledge of some of the orbital elements, for instance 
the spectroscopic ones (e, P and T). However, it appeared that 
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fixing u> at the value derived from the spectroscopic orbit of- 
ten led to orbital inclinations i unrealistically close to zero for 
spectroscopic binaries with radial velocity variations. Leaving 
the parameter u> free removes this difficulty, and offers more- 
over a way to check the consistency of the astrometric solution, 
since the astrometric co should be consistent with its spectro- 
scopic value. 

As far as outliers are concerned, we almost always keep 
the same dat a set as that used by FAST and/or NDAC (Vol. 
3, Sect. 17.6; [ESA (1997)| ). IAD that were not considered by 
FAST or NDAC (i.e., with the IA2 field set to T or 'n') were 
thus not included in our re-processing either. In a few cases, we 
noticed that because of the orbital contribution, some observa- 
tions yield residuals larger than 3<r of the residuals. In these 
cases, these observations were removed and the fit re-iterated. 
We never had to iterate more than twice to remove all outliers. 

2.1.2. Forcing parallaxes to be positive 

One of the five astrometric parameters entering x 2 m Eq. [I] is 
the para llax v j. With no other prescriptions as those described 
in Sect. |2.1.l| , the minimization process may very well end up 
with a negative parallax. Indeed, negative parallaxes are not 
rare in the Hipparcos and Tycho catalogues. 

Parallaxes cannot only be seen as inverse distances (which 
are defined positive) but also as the semi-major axis of the par- 
allactic ellipse (see Eq. fy. The direction of motion along the 
parallactic ellipse is of course imposed by the annual revolution 
of the Earth around the Sun, regardless of the actual dimension 
of the parallactic ellipse or of the observational uncertainties. In 
that sense, the parallactic ellipse is oriented, and negative paral- 
laxes can be seen as corresponding to a parallactic ellipse cov- 
ered in the wrong direction. That constraint being of physical 
nature, one should seek to fulfill it. In this section, we present a 
method which always delivers positive parallaxes. This method 
is especially useful for stars like those Mira variables or carbon 
stars that came out with large negative parallaxes in the Hippar- 
cos Catalogue. In those cases, forcing the parallax to be positive 
has a strong impact on the derived proper motion (Pourbaix et 
al., in preparation), which may be supposed to be better deter- 
mined with a physically- sound model yielding positive paral- 
laxes. However, the major drawback of the method is that the 
errors on the parallax do no longer follow a normal distribution. 
Therefore, the use of the parallaxes provided by this paper for 
e.g., luminosity calibrations should be done with care to avoid 
biases. 

In order to force the parallax to be positive - and at the same 
time avoiding the difficulties inherent to any constrained min- 
imization techniques - one may replace the constrained vari- 
able (the parallax vj) by an unconstrained oneQ. An appropriate 
choice appears to be 



which is C°° between ]0, +oo[ and M. The variable vj' is in 
fact equivalent to the distance modulus m — M since 



vj' = log w 



(7) 



M 



-5 — hvj . 



(8) 



The variable vj' is thus used in the minimization process in- 
stead of w as one of the p/. parameters entering Eq. [l] [with 
jtL, — jhL-zjln 10]. At the end of the minimization, Eq. [7] is 
reversed and the parallax w is derived from vj' . For any real 
number vj' , vj necessarily lies in between and oo. Such a 
substitution is legitimate, since the parallax is not a directly 
measured quantity, but rather one among many parameters used 
in a model fit to the observations. Moreover, since maximum- 
likelihood estimat orsPI enjoy the invariance property (see e.g., 
|Mood et al. 1974j p. 284), vj- and tu'-fitting must yield iden- 
tical results for those cases where tu-fitting yielded a positive 
parallax (since the vj' — log vj transform may be used when 
vj > to extract vj'). 

Most of the objects considered in the present paper have 
large parallaxes that would have come out positive by a direct 
fit of vj anyway. Thus, in the present case, the fitting of vj' 
(instead of m) does not represent so much of an improvement. 
Nevertheless, the procedure of ro'-fitting has been introduced 
here for the sake of generality. 

The price to pay is, however, that the errors on vj do not any 
longer follow a normal distribution. Moreover, the confidence 
interval of the parallax is no more symmetric. It has been esti- 
mated by the following expression: 



vj' — a < vj' < vj' + a 



(9) 



where vj' is the value resulting from the \ 2 - minimization, and 
a its estimated uncertainty. It should be noted that the errors on 
vj' do not follow a normal distribution either, since the normal- 



ity is only guaranteed for 'linear models' (Mood et al. 1974) 
and S (Eq. ^ does not depend linearly upon vj'. Therefore, the 
confidence interval corresponding to a given probability level 
is in general not symmetric around vj' . However, to ease the 
computations, the uncertainty a on vj' will be computed from 
the inverse of the Fisher information matrix at the point mini- 



mizing x ( see Sect. 2.1.3). This procedure implicitly assumes 



We are indebted to A. Albert for suggesting this trick 



that the model is linear in the vicinity of the minimum, so that 
the adopted confidence interval is in fact one that is symmetric 
around vj' . 

The second expression in Eq. || clearly shows that the par- 
allax is positive everywhere in the confidence interval, which 
would not be guaranteed in a constrained minimization of x 2 - 
That important property illustrates the superiority of this ap- 
proach with respect to the constrained minimization. 

2.1.3. x 2 minimization 

If the x 2 expressed by Eq. || were a quadratic expression of the 
unknown parameters pk, its unique minimum could be found 

3 x 2 is a maximum-likelihood estimator provided that the measure- 
ment errors are normally distributed 
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from the solution of a set of linear equations. However, the pa- 
rameters i, w and £1 enter x 2 m a highly non-linear way, so that 
the function expressing \ 2 m the 9-parameter space may have 
several local minima, and finding its global minimum is a much 
more arduous task. 

Faced with such situations, one of us (DP) has already suc- 
cessfully worked out global optimization techniques such as 
simulated annealing ( Pourbaix 1994 ; pourbaix 1998b| ). Prac- 
tical details about the implementation of the method to mini- 
mize the objective function \ 2 [Eq. p| | in the working space M 9 
may be found in Pourbaix (1998a)| . Simulated annealing be- 
ing a heuristic method, one can only prove its convergence to 
the global minimum after an infinite time (which we cannot af- 
ford). We thus stop the procedure after a finite time. In order 
to nevertheless have a good chance to obtain (a neighborhood 
of) the global minimum, we repeat 40 times this highly non 
deterministic minimization process. The best solution ever met 
(i.e., the one leading to the smallest x 2 value) is finally adopted. 
Once (a neighborhood of) the global minimum is thus obtained, 
it is tuned with the BFGS quasi-Newton algorithm (Dennis Jr. 
& Schnabel 1995). 

Unlike the Levenberg-Marquardt ( [vlarquardt 1963 ) mini- 
mization algorithm, BFGS does not return the covariance ma- 
trix of the model parameters. The inverse of the Fisher infor- 
mation matrix at the minimum is therefore used as the best es- 



timate of that covariance matrix ( [Pourbaix 1994| ). 

The whole procedure has been applied separately on the 
data from the FAST consortium only, from NDAC only and 
from both combined, thus resulting in three different solutions, 
hopefully consistent with each other. 

In a few instances, the solution obtained from the combined 
FAST+NDAC data set turns out to be very close to either the 
FAST or NDAC solution, but FAST and NDAC taken sepa- 
rately yield rather different solutions. That situation probably 
reflects the very different weights attributed to the two data sets 
for that particular object in the merging process applied to pro- 
duce the output catalogue. For our analysis we always keep the 
covariance matrices of the observations as they are given in the 
electronic version of the catalogue. 

As pointed out by an anonymous referee, in the case where 
all the Campbell elements (ao, i, lo and f2) are extracted from a 
fit to the astrometric data, they can advantageously be replaced 
by the Thiele-Innes elements (A, B, F and G) so that \ 2 be- 
comes a quadratic function of the model parameters. The min- 
imum of x 2 can then be found analytically and no minimiza- 
tion (neither global nor local) technique is needed. For the sake 
of generality, we nevertheless use the Campbell set (and thus 
the minimization scheme) because this more general scheme 
allows, if necessary, to easily incorporate external constraints 
(like for instance the knowledge of i for eclipsing binaries, or u> 
from the spectroscopic orbit ; see, however, the comment about 
fixing u> after Eq. |6| in Sect. [2.1. if ). Such additional constraints 
would be much more difficult to impose through the Thiele- 
Innes elements. 



2.1.4. F-test: 5-parameter vs orbital model 

The introduction of more free parameters in the orbital model 
as compared to the single star solution necessarily leads to a 
reduction of the objective function. To evaluate whether this 
reduction is statistically significant - or, equivalently, whether 
the orbital solution represents a significant improvement over 
the 5-parameter Hipparcos solution - requires the use of an F- 
test. 

The method used here is inspired from the test devised by 
Lucy & Sweeney (1971) If Xm P an< ^ Xorb denote the resid- 
uals for the Hipparcos 5-parameter model and for the orbital 
model (with 9 free parameters), respectively, the efficiency of 
the additional 4 parameters in reducing x^b below Xm p ma y 
be measured by the ratio: 



F = 



N-9 Xm P - xlrb 

Aorb 



(10) 



where N is the number of available measurements. 

If the hypothesis that there is no orbital motion (i.e., ao = 
in Eq. |5|) is correct, then it may be shown (Bevington & Robin- 
son 1992) that F follows a Snedecor F V1 U2 distribution with 
v\ = 4 and V2 = N — 9 degrees of freedom. Thus, if Eq. 
yields F = F, then, on the assumption that there is no orbital 
motion, the probability that F could have exceeded F is 



a = Prob(F > F) = Q{F\v u v 2 ). 



(11) 



In other words, a is the first risk error of rejecting the null hy- 
pothesis that the orbital and 5-parameter models are identical 
while it is actually true. 

The residuals given in the IAD files always relate to a 5- 
parameter solution, even when a more sophisticated model (the 
so-called 'acceleration' 7- or 9-parameter models, or even or- 
bital model) was published in the Hipparcos catalogue (see Ta- 
ble |lj). For those cases, the a value listed in Table [j] is always 
close to zero, although it does not really characterize the im- 
provement of the orbital solution with respect to the solution 
retained in the Hipparcos catalogue (which goes already be- 
yond a 5-parameter model). 

2.1.5. Astrophysical consistency of the orbital solution 

The minimization process will yield a solution in all cases, but 
that solution may not be astrophysically relevant. A statistical 
check of the significance of the orbital solution, based on the 
i^-test, has been presented in Sect. 2.1.4. In this section, two 



criteria testing the validity of the orbital solution on astrophys- 
ical grounds are presented. 

The first test is based upon the identity 



ao sin i P K 1 s/l — e 2 



2tt 



(12) 



where K\ is the radial-velocity semi-amplitude of the visible 
component. The left-hand side of the above identity entirely 
depends upon astrometric parameters, whereas its right-hand 
side contains only parameters derived spectroscopically. 
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This test has the advantage of being totally independent of tor grid: 
any assumptions. However, it involves the orbital inclination 
which is not always very accurately determined (see Sect. 3.2), 61 = 
so that the above identity may not always be very constraining 
considering the often large uncertainty on i. fr, 

A somewhat more constraining identity to assess the astro- 
physical plausibility of the computed astrometric orbit is the k 
following: 



Of) 

zu 



p2/3 



M 2 



(M x + M 2 ) 2 /3' 



(13) 



where M% and M 2 are the masses of the visible and invisible 
components, respectively. However, it relies upon assumptions 
regarding the stellar masses. In some cases (e.g., dwarf barium 
stars), the mass of the observable component may be estimated 
directly from the spectroscopically-derived gravity and from an 
estimate of the stellar radius from the spectral type (North et al. 
1999). For the other samples, an average mass derived from a 
statistical analysis of the spectroscopic mass function (Jorissen 
et al. 1998) is adopted. For all the samples considered here, 
the unseen component is almost certainly a white dwarf (the 
only possible exception being the Tc-poor S star HIP 99312 = 
HP 191589), whose mass may be taken as 0.62 ± 0.04 M 
dJorissen et al. 1998|). 



Because of the assumptions involved, this test is only used 
as a guide to identify astrophysically-unplausible solutions. It 
turns out that such cases are generally those with large error 
bars or with inconsistent N, F and A solutions, thus provid- 
ing further arguments not to retain those solutions. In very few 
cases (HIP 36042, 53763 and 60299), valid data yielded solu- 
tions not consistent with Eq. [□} Those cases were nevertheless 
kept in our final list. 



2.2. Transit Data 

Unlike the IAD, TD are only available for a small subset of the 
Hipparcos catalogue, e.g., for those stars that were known to be 
(or suspected of being) double or multiple systems at the time 
of the data reduction by the Hipparcos consortia. TD are a by- 
product of (or, more precisely, an input for) the multiple-star 
processing by the NDAC consortium. Another difference with 
respect to the IAD concerns photometry. Whereas IAD contain 
astrometric information only, the brightness of the (different) 
observable component(s) of the system can be retrieved from 
the TD. 



In the most general case (Quist & Lindegren 1999), each 
entry in the TD file corresponds to five numbers bi, . . . , 65 
which represent the coefficients of the first terms in the Fourier 
series modeling the observed signal as modulated by the detec- 



M Ii cos <f>i 

i 

-M^/jSin 4>i 

i 

Ar^J iC os(2^) 

i 

-iV^J iS in(2^ 



where Ii is the intensity of the i-th component of the system, 
and the phase (pi corresponds to its abscissa along the reference 
great circle. M = 0.7100 and N = 0.2485 are the adopted 
reference values for the modulation coefficients of the first and 



second harmonics (Vol. 1, Sect. 2.9; ESA 1997). In terms of the 
Cartesian coordinates (£,i,r/i) (see Eq. |3|) of component i in the 
plane tangent to the sky at the reference point specified in the 
TD file, the phase <pi writes 



4>i = Ixit + fyVi + fp("& ~ Wref)- 



(14) 



The reference point is assigned an arbitrary parallax nj rc f and 
proper motion, as given in the TD file. The phase derivatives 
f x , f y and f p with respect to £, r\ and w are also provided by 
the TD file. 

For the SB 1 systems we are interested in, the situation sim- 
plifies a lot since I2 may be taken equal to 0. The above system 
of equations is rank-deficient. From the second and third equa- 
tions, one can rewrite: 



4>= arg(b 2 , -b 3 ) 
and also express the uncertainty on <p as 



(15) 




'b 3 - 



With the above definitions, the objective function whose 
minimum is sought is given by: 



1 



N 



D =N^ 

k=l 



\4>k -^fe) 2 



(16) 



where N is the number of TD. In the above expression, 4>k is 
the observed phase at time tk as derived from Eq. |5l whereas 
<pk is the phase computed from Eqs. || and [l4] (noting that 
f x RP a + fyRPs = f P ) for a given set of astrometric and or- 
bital parameters. 

The remaining of the method follows the same steps as de- 
scribed in relation with the IAD, i.e. global and local optimiza- 
tion, positiveness of the parallax, . . . 
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3. Results 



Table [T] lists the various samples that have been considered, 



namely dwarf barium and subgiant CH stars (McClure 1997 



North, priv. comm.), mild and strong barium stars (Jorissen 
et al. 1998), Tc-poor S stars (Jorissen et al. 1998) and CH stars 
(McClure & Woods worth 1990). Their spectroscopic orbital 
parameters have been taken from the reference quoted. 

Table [j] also provides various parameters that allow either 
to assess the quality of the derived orbital solution or to identify 
the reason why a reliable orbital solution could not be derived 
for some of the spectroscopic binaries considered. The follow- 
ing parameters potentially control our ability to derive an astro- 
metric orbit from the IAD, with favorable circumstances being 
mentioned between the parentheses: the parallax (column 3; 
large parallax), the number of available IAD measurements 
(FAST+NDAC; column 4; large number of measurements), the 
ecliptic latitude (column 5; this parameter may play a role since 
it controls how different the orientations of the reference great 
circles are; favorable cases have absolute values larger than 
47°), the orbital period (column 6; in the range 1 - 3 yr to 
ensure a good sampling of the orbit), the orbital eccentricity 
(column 7; low eccentricity). 

The following columns characterize the quality of the solu- 
tion obtained from the minimization process. Column 1 1 pro- 
vides x 2 / {Niad — 9), which should be of the order of unity if 
the internal error a (Eq. |^) on the abscissae has been correctly 
evaluated by the reduction consortia and if the model provides 
an adequate representation of the data. The first risk error as- 
sociated with rejecting the null hypothesis that the orbital and 
Hipparcos 5-parameter solutions are identical (Eq. [IT]) is given 
in column 12. Low values of a are generally associated with 
Hipparcos solutions of the G, X or O types (as listed in col- 
umn 9; see the caption to Table [j] for more details), since the 
orbital motion is then large enough to have been noticed al- 
ready by NDAC or FAST. Columns 13 and 14 compare the 
aa/zu ratio derived from the orbital solution to its expected 
value from Eq. [l3|. In columns 11 to 13, the data refer to the 
orbital solution obtained by combining NDAC and FAST data. 
For dwarf barium stars, the masses used to estimate ao jw ac- 
cording to Eq. H are M 2 = 0.62 ± 0.04 M , whereas M x is 
derived from the spectroscopic gravity, with an estimated er- 
ror of 0.05M Q (North, priv. comm.). According to the statisti- 
cal analysis of the spectroscopic mass functions performed by 
Jorissen et al. (1998), Mi and pairs (expressed in M Q ) of 



(1.7 ± 0.2, 0.62 ± 0.04), (2.1 ± 0.2, 0.62 ± 0.04) and (1.8 ± 
0.2, 0.62 ± 0.04) have been adopted for strong barium stars, 
mild barium stars and Tc-poor S stars, respectively. The same 
analysis performed by McClure & Woods worth (1990)| for CH 
stars yielded Mi = 1.0 ± 0.1 M and M 2 = 0.62 ± 0.04 M . 

Astrometric orbits were accepted when a is smaller than 10 
per cent and the expected a^jw value falls within the la con- 
fidence interval. A few cases not fulfilling these criteria were 
nevertheless accepted after visual inspection of the orbital arc. 

Examination of Table [j] reveals that the following criteria 
need to be fulfilled in order to be able to extract a reliable astro- 



metric orbit from the Hipparcos data: w <; 5 mas, P 10 yr, 
N <; 25. The success rate is as follows for the various samples: 
dwarf barium stars (9/13), mild barium stars (3/24), strong bar- 
ium stars (6/26), Tc-poor S stars (3/10) and CH stars (2/8). The 
high success rate for dwarf barium stars naturally results from 
the fact that these dwarf stars are on average closer from the 
sun than the giant stars. 

Table ^| lists the astrometric and orbital parameters for the 
reliable orbits according to the criteria discussed above. The 
results from the different processing modes are collected in Ta- 
ble ||, according to the symbol given in column 2: H/S refers 
to the parameters from the Hipparcos catalogue and from the 
spectroscopic orbit (on that line a,Q /w is the semi-major axis in 
A.U. estimated from Eq. [jj]), F refers to the processing of the 
IAD from FAST only, N to the IAD from NDAC only, A from 
the processing of the combined FAST/NDAC data set, O to the 
orbital parameters from the DMSA/O, and T to the parameters 
resulting from the TD. 

Most of the retained orbits are indeed characterized by 
X 2 / (N — 9) values of the order of unity, as expected. The first- 
risk errors a are not always close to 0, but if the derived value 
for the semi-major axis is in good agreement with its expected 
value, that agreement has been considered as sufficient for re- 
taining the orbit. The only cases where the reverse situation 
occurs (small a but discrepant ao/xu) are the dwarf barium star 
HIP 60299, the mild barium star HIP 36042 and the CH star 
HIP 53763. Although the orbit of the latter is not well defined, 
it has been kept in our final list to illustrate the large uncertainty 
on w resulting from an orbital period close to 1 yr (Sect. 

The F, N and A solutions for the retained orbits are also 
generally in good agreement, the only exceptions being the 
dwarf Ba stars HIP 62409 and HIP 116233, and the mild Ba 
star HIP 1 17607. However, the model parameters of these sys- 
tems are highly correlated, and the different measurement er- 
rors in the different data sets thus drive the solution in different 
directions. This statement may be expressed in a quantitative 
way using the concept of efficiency e in troduced by Eichhorn 
(1989) and |Pourbaix & Eichhorn (1999)|. It is defined as 



n 



k=l 



rife=i ikk 



(17) 



where are the eigenvalues of the covariance matrix of the 
estimated parameters, qkk are its diagonal elements, and p de- 
notes the number of parameters in the model. If e is close to 
unity, there is obviously little correlation between the param- 
eters. For the combined NDAC+FAST solution, it amounts to 
0.32, 0.21 and 0.42 in the case of HIP 62409, HIP 116233 and 
HIP 1 17607, respectively, thus translating some degree of cor- 
relation between the model parameters. 

In Table 0, the uncertainty on ao jvj has been computed by 
combining the upper and lower limits on ao and m, thus ne- 
glecting any possible correlation between these two quantities 
(which is generally small - except for the three systems listed 
above - as derived from the efficiency being close to unity). 

The orbital solutions derived in the present paper are too 
many to display the astrometric orbit for all cases. A few repre- 
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Table 1. The various samples of CPRS that have been considered in the present work, along with various parameters charac- 
terizing the quality of the orbital solution derived from the IAD (see text for details), tuhip is the parallax from the Hipparcos 
catalogue, Njad is the number of IAD available for the star considered, j3 is the ecliptic latitude, T is the time of periastron 
passage (or of maximum velocity for circular orbits). SoIhip is the number of parameters used in the Hipparcos solution. It may 
be one of the following types: '5' = 5-parameter solution, '7' = 7-parameter solution, i.e., including a linear acceleration term 
in the proper motion (see Sect. 2.3, Vol. 1 of ESA 1997), '9' = 9-parameter solution, i.e., including a quadratic acceleration 



term in the proper motion ('7'- and '9'-entries are listed in the 'G' section of the Double and Multiple Star Annex - DMSA - 
of the Hipparcos and Tycho Catalogues), 'X' = stochastic solution (described in DMSA/X), and 'O' = orbital solution given in 
the DMSA/O. The column labeled 'Transit?' indicates whether TD are also available for the corresponding star. The column 
labeled x 2 /(-^iad — 9) lists the unit- weight variance for the orbital-parameter model, a is the first risk error of rejecting the null 
hypothesis that the orbital and single star models are identical while it is actually true (Eq. [n]). The expected value of the orbital 
separation ao/vj is computed from Eq. (|l3|). The lcr confidence interval is given next to the computed astrometric ao/m value. 
Column 'Rem' indicates whether the orbital solution derived from the IAD has been accepted. An asterisk in column 'Rem' 
refers to a note at the bottom of the Table 
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Remarks: 



HIP 53763 & 62827: large uncertainty on vo in the 9-parameter solution (see Table 0), since the orbital period is close to 1 yr. 



sentative cases among the different subsets of Table (orbital 
periods shorter or longer than the duration of the Hipparcos 
mission, small or large parallaxes . . . ) have instead been se- 
lected and are presented in Fig. [j]. 



3.1. Comparison with DMSA/O solutions 

For HIP 17296 (Tc-poor S), 31205 (strong Ba) and 56731 
(strong Ba), orbital solutions are provided in the DMSA/O and 
were derived using spectroscopic elements from the literature. 
HIP 50805 (dwarf Ba) is the only case in our sample where an 
orbit could be derived from scratch by the Hipparcos consor- 
tia. For all these systems, the astrometric orbits derived by the 
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Table 1. (Continued.) 
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Rem: HIP 1 12821 is also sometimes classified as a dwarf barium star 



methods described in Sect. || are in excellent agreement with 
the DMSA/O elements, thus providing an independent check 
of the validity of our procedures. Further checks are presented 



in Sect. 3.2. 



The large number (23) of systems for which orbital solu- 
tions could be extracted from the Hipparcos data (as compared 
to only 4 of those already present in the DMSA/O) illustrates 
the great potential that still resides in the Hipparcos IAD or TD. 

3.2. Check of the astrometric orbit 

Several checks are possible to evaluate the accuracy of the as- 
trometric elements derived in the present paper. 

First, it is possible to compare the astrometric and spectro- 
scopic values of to, the argument of periastron. In most cases, 
the two determinations agree within 2a (Fig. ||). However, even 
when the orbital period is shorter than the Hipparcos mission, 



the lo derived from the IAD is seldom as precise as the spectro- 
scopic one. 

In a few cases, the spectroscopic orbit is assumed to be cir- 
cular. In that case, the time T of passage at periastron becomes 
meaningless, and is replaced by the time of the nodal passage 
(or, equivalently, the time of maximum radial velocity). This is 
equivalent to setting lo equal to 0. Non-zero values for lo would 
correspond to other conventions for the origin epoch. Among 
our systems with acceptable orbits, two have circular orbits: 
HIP 53763 (CH star) and HIP 52271 (strong barium star). For 
these systems, our fit leads to values of lo significantly different 
(2a) from 0, indicating that at the reference epoch, the star is 
in fact far from the node where it was expected to be. 

It is also possible to compare the astrometric value of K\ 
(using in Eq. [l~2| the value of a, i and w from the astrome- 
try and e and P from the spectroscopy) with the spectroscopic 
value. Fig. [3] shows that, even if a^jw is well defined, the in- 
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1.17±0.42 


Accepted 


15264 


20394 


2.16±1.14 


48 


-15.3 


2226.0 


0.2 


47929 


5 




1 .39 


77 


8 


1.25±0.44 




1 7402 


24035 


3.72±0.8 


70 


-76.9 


377.8 


0.02 


48842.65 


5 




1.12 


49 


3.3e+09 + °° 


0.38±0.14 


* 


23 16H 


3 1487 


4.54±1.21 


40 


294 


1066.4 


0.045 


45173 


5 




0.45 


86 


0.47 J' 7 


0.77±0.27 




25452 


36598 


3.32±0.68 


72 


-85.3 


2652.8 


0.084 


45838.95 


5 




0.99 


5 


« 11 


1.41±0.50 




zyoyy 


42537 


-1.13±0.93 


64 


-75.9 


3216.2 


0.156 


46147.32 


5 


y 


1 .30 


82 


1.3e+15 + °° 


1.60±0.57 




29740 


43389 


-1.25±1 


52 


-25.8 


1689.0 


0.082 


47222.46 


5 




0.53 





^ 2400 
" 0.4 


1.04±0.37 


* 


30338 
31205 


44896 
46407 


1.56±0.7 
8.25±0.92 


90 
80 


-56.9 
-34.3 


628.9 
457.4 


0.025 
0.013 


48464.3 
47677.45 


5 
O 


y 


0.88 
0.77 







1 8 

0.92 
°- 71 0.46 


0.54±0.19 
0.44±0.15 


Accepted 


32713 


4964 1 


0.73±0.88 


48 


-19.2 


1768.0 





46306 


5 




0.63 





7 5 22 
3 2.3 


1.08±0.38 




32960 

"iftft/Ti 

50006 


50082 
ou i y / 
88562 


4.71±0.99 
1.67±0.84 
3.13±1.17 


40 
92 
50 


-16.2 
-50.5 
-25.1 


2896.0 
3243.8 
1445.1 


0.188 

0.34 

0.204 


45953.12 
46015.97 
45781.71 


5 
5 
5 




1 .39 
0.64 
0.95 


44 
59 
4 


, „ 8.8 
J - y 0.62 

21 69 
66 

6.5 

i - L 1.3 


1.49±0.53 
1.61±0.57 
0.94±0.33 




5227 1 


92626 


3.4±0.71 


104 


-50.5 


918.2 





49147.83 


5 




1 .05 


9 


99 
0.59 £jg 


0.69±0.25 


Accepted 


56404 
56731 


100503 
101013 


0.67±1.2 
7.07±0.68 


64 
74 


-30.8 
43.2 


554.4 
1710.9 


0.061 
0.195 


46144.83 
43934 


5 
O 




0.63 
1.19 


26 



23 21 
15 

1.4 

0.88 


0.50±0.18 
1.05±0.37 


Accepted 


60292 


107541 


5.78±1.36 


68 


-29.5 


3569.9 


0.104 


44388.16 


5 




0.96 


2 


27 51 
11 12 


1.72±0.61 




68023 


121447 


2.21±1.02 


52 


-6.0 


185.7 





46922.35 


5 


y 


1.14 


19 


3 5 7 6 
J - 3 1.3 


0.24±0.08 




69290 


123949 


.97±1.32 


34 


-5.6 


9200.0 


0.972 


49144.96 


5 




1 .50 


92 


470 9000 


3.23±1.14 




94103 


178717 


2.9±0.95 


52 


32.5 


2866.0 


0.434 


44258 


5 


y 


0.89 


71 


16 Q 


1.48±0.53 




101887 


196445 


1.49±1.54 


50 


-21.2 


3221.3 


0.237 


46037.95 


5 


y 


1 .06 


63 


9.3 I 00 


1.60±0.57 




103546 


199939 


3.16±0.75 


72 


57.6 


584.9 


0.284 


45255.1 


5 




0.79 


26 


0-44 ° H 


0.51±0.18 


Accepted 


104542 


201657 


4.49±1.07 


64 


31.6 


1710.4 


0.171 


46154.95 


5 




1 .06 


73 


9 

0.81 q 


1.05±0.37 




104684 


201824 


0.56±1.56 


50 


7.4 


2837.0 


0.342 


47413 


5 




1.79 


87 


2.1e+03 + °° 


1.47±0.52 




110108 


211594 


4.59±1.18 


26 


4.4 


1018.9 


0.058 


48538.19 


5 




1.11 


17 


1 Q 

°" a43 


0.74±0.26 


Accepted 


Tc-poor S stars 




























5772 


7351 


3.21±0.82 


28 


19.1 


4592.7 


0.17 


44696 


5 


y 


0.92 


88 


ft a 57 
6.9 


1.97±0.67 




8876 


+21°255 


-1.92±1.5 


38 


9.5 


4137.2 


0.209 


43578.31 


5 


y 


1.13 


87 


2.4e+13 + °° 


1.84±0.66 




17296 


22649 


6.27±0.63 


78 


42.2 


596.2 


0.088 


42794.5 


O 


y 


1.22 


1 


0.52 
0.3/ q 24 


0.51±0.17 


Accepted 


25092 


35155 


1.32±0.99 


52 


-31.7 


640.6 


0.071 


48092.41 


5 


y 


0.79 


51 


°- 94 lis 


0.53±0.18 


Accepted 


32627 
38217 


49368 
63733 


1.65±1.11 
0±0.99 


60 
84 


-17.4 
-39.2 


2995.9 
1160.7 


0.357 
0.231 


45145.37 
45990.92 


5 
5 


y 


1 .03 
0.95 


57 
58 


1ft 110 

16 0.91 
1.76+14+°° 


1.48±0.50 
0.79±0.27 




90723 


170970 


1.83±0.67 


86 


59.4 


4392.0 


0.084 


48213.21 


5 


y 


0.98 


39 


27 59 
11 8.8 


1.91±0.66 




99124 
99312 


191226 
191589 


0.39±0.71 
2.25 ±0.77 


84 
84 


54.9 
52.0 


1210.4 
377.3 


0.19 
0.253 


49691.78 
48844.02 


5 
5 




0.98 
0.80 


23 
82 


11 620 
11 0.14 

2 4 
0.58 q 


0.81±0.27 
0.37±0.13 


Accepted 


115965 


+28°4592 


1.72±1.26 


64 


29.6 


1252.9 


0.091 


48161.32 


5 




0.78 


72 


i.4 r 


0.83±0.28 





HIP 29740: a = but the solutions are totally unrealistic 

HIP 17402: vj not accurately determined, since orbital period is close to 1 yr. 



clinations are generally not very accurately determined, thus 
leading to uncertain values of K\ . This unfortunate property of 
i is well illustrated on Figs. |] and ^. One example where the 
accuracy of the astrometric value of i must be questioned is the 
dwarf barium star HIP 105969: despite the fact that the astro- 



metric ao/zu ratio perfectly agrees with its estimate based on 
the masses, the astrometric prediction of the semi-amplitude of 
the radial-velocity variations differs by almost of factor of 2 as 
compared to the actual spectroscopic value (Table ||). The only 



Table 2. Astrometric orbital elements for various classes of CPRS. The first line provides the data from the HIPPARCOS catalogue or from the spectroscopic orbit (label 
H/S in column 2). The following lines provide the orbital solutions derived in the present work from the FAST, NDAC and combined FAST+NDAC IAD (labels F, N and 
A, respectively in column 2). When available, the orbital solution from the transit data (label T in column 2) or from the DMSA/O (label O) are given next. The columns 
zu, a /zu and K\ contain the nominal value followed by the boundaries of the Tcr' confidence interval (computed by propagating the la errors on a , zu' and i as if they 
were uncorrected). An asterisk after the HIP number refers to a note at the end of the Table 



HIP 



Sol 



zu 
(mas) 



Ma 

(mas/y) 



(mas/y) 



do 
(mas) 



i 

o 



o 



O 



ao/zu 
(AU) 



(km/s) 



N 



X X a 
HIP Orb (%) 



CH stars 

53763* H/S 
F 

N 

A 

104486 H/S 
F 



5.29±1.47 

o oo 8-51 
2.88 o 97 

, 5 o 15.37 
LM 0.16 
,8.18 



2.70 



0.89 



N 
A 

Dwarf Ba stars 

32894 H/S 
F 

N 

A 



4.93±0.84 
5.63 
3.83 
7.67 
5.60 
6.01 
4.32 



4.64 
6.56 
5.09 



14.11±1.96 
15.41 
13.04 
15.03 
12.24 
15.45 
13.15 



14.17 
13.57 
14.25 



49166 



H/S 
F 

N 

A 



7.9±1.39 
;.62 



6.94 
7.40 
7.26 



5.59 
9.27 
5.91 
8.83 
5.97 



-24.16±1.18 
-23.28±1.3 

-24.22±1.6 

-23.61±1.2 

-43±0.66 
-48.91±2.9 

-48.78±3.2 

-48.28±2.7 



34.81±1.0 
35.16±1.2 
35.07±0.9 

-65.45±0.85 
-66.42±1.5 

-66.24±1.6 

-66.38±1.4 



-15.03±1.2 
-14.66±1.2 

-15.43±1.5 

-15.13±1.2 

29.96±0.52 
32.21±2.1 

32.23±2.4 

32.43±2.0 



32.62±1.36 106.64±1. 



108.50±1.4 
108.49±1.8 



4.92±1.9 
5.21±2.5 
4.91±1.9 

8.18±1.4 
7.34±1.5 
7.46±1.2 



12.29±1.2 
11.43±1.4 



108.41±1.4 11.99±1.2 
70.1±1.1 

70.59±1.2 3.81±1.6 

70.45±1.5 4.42±1.8 

70.50±1.2 3.44±1.5 



62.6±22.9 
79.4±19.7 
68.7±19.3 

71.3±19.0 
66.7±24.9 
64.8±21.6 



108.4±6.5 
111.7±8.7 
109.0±6.5 



113.0±44.3 247.6±44.1 1.71 
100.2±34.7 

103.5±37.6 254.2±36.0 



0.44±0.22 

6.98 



0.36 

258.4±32.6 3.27 Q^g 3 
7.67 



1.81 



0.36 



133±20 
74.7±22.2 

86.5±30.6 

75.7±24.5 

223±12 
226.8±6.5 

219.8±8.2 

225.1±6.2 



128±3 
113.6±33.1 107.7±30.7 

88.8±28.1 87.8±23.9 

109.1±32.3 99.7±30.6 



103.1±10.6 
97.7±16.3 
104.4±13.6 



69.7±8.7 
64.1±12.1 
69.2±8.9 

341.6±23.2 
356.7±15.1 
338.1±21.3 



1.2±0.58 
2.49 



1.76 
1.12 
1.46 



1.21 
1.58 
0.76 
2.01 
1.04 



0.89±0.46 
1.03 



0.87 
0.84 
0.84 



0.72 
1.05 
0.67 
1.00 
0.70 



0.39±0.19 
0.97 



0.55 
0.60 
0.47 



0.26 
1.04 
0.29 
0.83 
0.22 



5.42±0.28 
230.64 



50.25 
106.60 
56.00 



7.62 

1565.08 

5.08 

254.06 

9.16 



5.3±0.19 
18.62 
7.15 
11.76 
3.80 
14.98 
5.30 



12.46 
7.67 
9.89 



9.7±0.2 
12.11 



9.85 
9.38 
9.53 



7.83 

12.27 

6.87 

11.70 

7.58 



12.28±0.1 

on 1 c 39.08 
20.36 5 73 

42.24 



24.15 
18.11 



10.16 
33.52 
5.54 



67 



32 


33.7 


22.8 





33 


28.9 


21.0 





70 


61.2 


51.4 


1 


99 








48 


90.2 


69.2 





49 


68.4 


49.1 





100 


145.0 


123.3 





77 








38 


159.4 


37.6 





39 


113.8 


40.8 





78 


191.3 


72.2 





87 








42 


46.5 


40.5 


3 


43 


48.8 


42.2 


2 


90 


88.2 


83.4 


34 



Table 2. (Continued.) 



HIP 



Sol 



(mas) 



(mas/y) 



(mas/y) 



a 
(mas) 



(°) 



(AU) 



(km/s) 



N 



X 

HIP 



X 

Orb 



t/2 
§ 

u 

CO 



I 

o 

I 

S 

■c 
o 

< 

X 
•5 

■e 
§ 

Oh 



Dwarf Ba stars (cont.) 

50805 H/S 
F 

N 



60299 



62409* 



69176 



A 
O 
T 

H/S 
F 

N 

A 

H/S 
[F] 
N 
[A] 

H/S 
F 

N 

A 



23.42±0.93 
23.43 24 - 36 



22.94 
23.40 



23.5 



22.52 
24.01 
21.91 
24.28 
22.54 

24.6 
22.3 



5.02±1.06 
- 6.59 



5.37 
6.02 
5.64 



4.37 
7.38 
4.91 
6.77 
4.69 



8.2±1.28 
9.31 
6.30 
10.18 
6.56 
9.25 
6.38 



7.66 
8.17 
7.68 



8.75±1.39 
12.46 



10.36 
9.80 
10.16 



.61 
12.31 
7.79 
12.19 
8.46 



-26.77±0.67 
-27.78±0.8 

-25.91±0.8 

-26.89±0.7 

-26.3±0.83 

-147.79±0.96 
-141.76±2.1 

-139.47±2.6 

-141.41±1.9 

-136.5±1.63 
-137.04±7.9 

-126.05±8.9 

-136.08±7.2 

34.08±1.22 
35.36±1.6 

33.08±1.8 

34.38±1.5 



55.41±1.36 
56.42±0.9 

56.12±1.1 

56.51±0.9 

56.4±0.98 

9.32±0.7 
10.25±1.9 

11.74±1.5 

11.20±1.2 

-19.28±0.79 
-15.07±3.7 

-16.19±3.8 

-15.09±3.3 

-29.9±1.12 
-28.80±1.2 

-30.76±1.4 

-29.26±1.1 



15.75±1.0 

15.28±1.2 

15.59±0.9 
15.76±1.0 
16±1.3 

9.35±2.5 
12.02±3.0 
9.95±2.4 

5.94±4.1 

11.88±8.8 

6.71±4.7 

4.64±1.1 
5.07±1.4 
4.81±1.1 



102.1±3.0 

99.8±3.3 

100.6±2.8 
100.2±2.9 
99±3.6 

116.9±13.1 
118.0±10.0 
123.3±11.6 

81.4±42.6 
120.8±32.5 
94.2±34.1 

52.6±31.0 
63.8±27.6 
53.9±28.4 



137±4 
144.7±4.5 

147.7±5.5 

146.3±4.3 



142±5.3 




69±48 




127.5±13. 


.3 


113.9±13. 


.0 


124.1±14. 


.3 


186±24 




167.3±25. 


.3 


276.0±35. 


.5 


168.2±21. 





192±18 




203.8±32. 


.2 


191.5±22. 


.9 


203.5±28. 


.2 



343.1±4.0 

343.2±4.4 

343.8±3.7 
344.1±3.6 
342±4.6 

266.6±18.9 
258.5±15.0 
264.1±16.6 

182.8±80.2 
293.1±31.2 
165.1±57.7 

17.0±40.9 
30.9±30.2 
8.5±35.5 



0.67±0.3 
0.74 



0.67 
0.67 

°- 67 0.60 
0.67±0.05 
0.78 



0.61 
0.75 
0.59 
0.73 



0.69 



0.61 



0.83±0.26 
2.72 



1.74 
1.99 
1.76 



1.04 
3.05 
1.22 
2.63 
1.12 



9.67 ±0.1 
12.04 



10.79 
10.77 
10.75 



9.60 

12.23 

9.39 

11.92 

9.63 



2.32±0.16 
21.33 



12.56 
14.24 
11.93 



6.41 

23.47 

7.80 

19.74 

6.40 



1.3±0.65 7.17±0.22 



0.78 
1.45 
0.87 



1.60 
0.19 
3.15 
0.30 
1.79 
0.22 



0.44±0.16 
0.67 



0.45 
0.52 
0.47 



0.28 
0.84 
0.29 
0.69 
0.31 



4.69 
7.64 
5.33 



9.80 
0.74 
19.26 
0.84 
10.95 
1.04 



12±0.43 
8.48 * 5 J 3 

11.05 \ g if 
16.37 



58 



9.11 



3.15 



29 


296.4 


22.2 





29 


202.9 


16.3 





60 


330.2 


42.0 





98 








52 








25 


62.4 


15.3 





25 


57.6 


15.9 





50 


96.3 


38.1 





57 








28 


21.9 


15.6 





29 


36.3 


17.5 





58 


70.2 


63.0 


25 


55 








25 


46.1 


23.6 





27 


37.3 


15.6 





54 


82.3 


55.9 






Table 2. (Continued.) 



HIP 



Sol 



w 
(mas) 



Ma 

(mas/y) 



H6 
(mas/y) 



fflo 
(mas) 



i 

(°) 



(°) 



(°) (AU) (km/s) HIP Orb (%) 



Dwarf Ba stars (cont.) 



105969 



107818 



116233* 



118266* 



H/S 
F 

N 

A 

H/S 
F 

N 

A 

H/S 
F 

N 

A 

H/S 
F 

N 

A 



16.61±1.78 



16.10 
15.92 
16.22 



16.96 
15.29 
16.84 
15.05 
17.01 
15.46 



7.53±1.51 
10.54 
7.12 
10.79 
6.92 
10.61 
7.30 



8.66 
8.64 
8.80 



8.83±1.21 
11.29 
8.43 
12.34 
8.56 
10.97 
8.35 



9.76 ; 



10.28 
9.57 



6.95±1.45 
10.90 
6.72 
10.35 
5.35 
10.84 
6.76 



8.56 
7.45 
8.56 



200.64±1.56 
200.75±0.8 

200.27±0.7 

200.25±0.6 

16.16±1.24 
16.50±1.6 

16.89±1.8 

16.95±1.5 

27.08±1.69 
26.61±4.0 

44.68±17.8 

26.58±3.8 

-30.31±1.37 
-29.70±1.5 

-30.71±1.7 

-30.08±1.4 



-93.14±1.52 
-93.48±0.7 

-92.23±0.8 

-93.07±0.6 

-36.55±0.8 
-36.90±1.0 

-36.43±1.1 

-36.35±1.0 

-22.19±1.14 
-9.47±3.3 

-16.68±6.3 

-9.98±3.2 

-121.84±0.86 
-120.64±1.1 

-121.93±1.3 

-120.97±1.0 



12.27±0.9 132.7±6.9 

12.07±1.0 132.1±7.5 

12.12±0.8 133.1±6.6 

3.09±1.4 107.0±40.5 

4.23±2.3 126.7±31.0 

3.44±1.4 107.3±32.4 

11.10±1.4 26.7±50.4 

14.91±9.0 76.8±17.4 

11.11±1.4 23.0±54.2 

2.60±1.7 45.2±66.7 

3.11±2.1 112.1±29.6 

2.30±1.4 60.1±55.6 



192±21 
225.9±11.2 

239.6±12.0 

233.0±10.8 

281 ±42 
306.0±27.3 

340.0±55.9 

302.1±28.5 

196±17 
355.1±91.1 

312.1±36.2 

201.9±118.9 

250±44 
360.6±73.3 

267.6±48.7 

205.2±47.3 



273.6±10.9 
282.7±11.8 
277.5±10.6 

153.4±43.2 
226.8±48.4 
175.7±37.3 

158.1±87.7 
249.5±32.7 
312.0±115.1 

219.1±68.2 
57.0±38.7 
18.5±53.7 



0.77±0.33 



0.76 



0.86 



1 0.67 

°- 76 0.66 
0.84 



0.75 



0.66 



0.58±0.22 
0.63 
0.16 
0.95 
0.18 
0.66 
0.19 



0.36 
0.49 
0.39 



0.97±0.37 
1.48 
0.86 
2.79 



1.14, 



1.45 
1.16 



0.48 
1.50 
0.89 



0.4±0.19 
,0.64 
0.09 
0.98 
0.10 
0.55 
0.08 



0.30 
0.42 
0.27 



3.29±0.15 

700 1:11 

,90 
5.34 
8.41 



7.03 



6.81 



5.36 



10.4±0.19 
10.22 
1.41 
15.29 
1.09 
10.71 
2.05 



5.54 
6.37 
6.06 



7.42±0.26 
11.24 
0.00 
,21.75 



3.98 , 



11.00 
3.53 



3.22 
11.38 
0.00 



8.55±0.39 
c/x 22.50 

'■ bi 0.00 

B n i 34 . 53 
■ n 2.12 
19.46 



8.24 



0.23 



77 



38 


334.0 


24.7 





39 


290.9 


40.2 





80 


418.7 


72.1 





50 








25 


22.5 


16.3 


1 


25 


24.8 


15.4 





52 


42.7 


35.1 


7 


54 








23 


149.7 


4.6 





27 


113.4 


17.6 





56 


209.5 


39.3 





79 








39 


42.2 


37.3 


7 


39 


25.9 


22.6 


5 


80 


74.7 


71.0 


45 



Remarks: 

HIP 62409: Only NDAC yields a realistic orbit. 
HIP 116233, HIP 118266: Orbits not well denned. 

HIP 53763: orbit not well denned with large uncertainty on zu since P is close to 1 yr. 
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Table 2. (Continued.) 

HIP Sol w /4 T0& Qfl i w ^ ao/zu i?i N a 



(mas) (mas/y) (mas/y) (mas) Q Q Q (AU) (km/s) HIP Orb (%) 



Mild Ba stars 




























36042* 


H/S 
F 

N 

A 


2.36±0.97 

ZJii 1.89 

2 00 L13 
2.44 3|9 


2.74±0.94 
2.49±1.0 

2.78±1.1 

2.60±1.0 


-0.74±0.61 
-0.51 ±0.7 

-0.20±0.7 

-0.31±0.6 


5.79±1.2 
3.57±1.3 
4.86±1.1 


94.0±9.7 
119.1±22.7 
101.0±11.5 


17±4.6 
81.8±13.3 

61.4±30.5 

75.2±14.9 


289.1±8.1 
296.7±23.4 
295.8±9.9 


0.5±0.15 
2.08 fl\ 

1 78 4:32 
lJ * 0.64 

1.99 HI 


6.91±0.14 
34.39 fl-53 

25-86 6 71 5 ^ 8 

qo nc 61.16 

32 - 36 15.56 


56 
28 

28 

56 


33.7 
28.0 
63.1 


15.8 
18.5 
46.9 




1 


105881* 


H/S 
F 


8.19±0.9 

9 05 10 - 23 
yu 5 8.00 


-2.61±0.96 
-0.55±6.6 


18.88±0.42 
26.59±3.4 


9.77±4.7 


102.8±17.9 


239.51±10 
228.6±12.5 


187.6±45.2 


1.2±0.35 
1-08 0J0 


2.763±0.161 

5 0? 8 - 62 
5.02 2 04 


51 

25 


29.8 


18.9 







N 


7 71 9.08 
6.54 


4.71±7.4 


23.37±3.8 


9.24±8.3 


119.8±31.2 


269.4±32.3 


237.5±45.1 


1-20 §:?f 


4.96 2 2 F 


26 


24.2 


12.8 







A 
T 


a fin 9-67 
8 - 60 7.65 

7 7 7.8 

'•' 7.5 


-0.57±5.8 
9.5±0.86 


25.30±3.1 
24.5±0.49 


8.56±4.3 
12.3±0.89 


111.7±20.7 
112±2.6 


233.8±19.8 
261±2.7 


190.7±50.0 
240±3.5 


1 fi 1-75 
16 1.46 


4 - 4 i 


54 
86 


50.2 


41.4 


6 


117607* 


H/S 
F 


4.61±0.95 
4« 5.59 
4 - 35 3.39 


2.88±1.29 
6.43±3.6 


3.32±0.92 
0.57±1.4 


3.55±1.7 


116.0±38.8 


157.61±25 
139.1±56.0 


128.7±57.5 


0.78±0.23 

0-82 o:i 


3.63±0.08 

, ,0 13.13 
6.18 L17 


60 
29 


104.1 


18.5 







N 
A 


3 q i 5.40 
j.yi 2.84 

t Q7 5.07 
J - y ^ 3.03 


-0.52±3.2 
3.33±2.7 


-3.07±1.6 
-0.58±1.3 


6.83±2.6 
3.41±1.0 


83.8±15.5 
82.3±37.1 


138.4±21.1 
154.6±23.4 


214.8±18.6 
173.3±42.4 


!- 74 0:79 

87 1M 
U.8/ Q 47 


14 - 61 6.18 
7 77 12.30 
2.83 


31 

62 


111.5 
150.2 


15.0 
39.9 







Remark: 



HIP 36042: large discrepancy between astrometric and astrophysical ao jvj. However, the orbit has been retained since the correlation between the 9 parameters is very low (the average 

efficiency of the three solutions is 0.69; see Eq. 17). 
HIP 105881: orbit not well defined. 

HIP 1 17607: although the orbit is not well defined, it has been kept in the Table to illustrate the correlations that may appear between ao and /j, 



Table 2. (Continued.) 



HIP 



Sol 



(mas) 



Ha 

(mas/y) 



H6 

(mas/y) 



a 
(mas) 



i 

n 



o 



H K[ N x 2 x 2 ~ 
(°) (AU) (km/s) HIP Orb (%) 



Strong Ba stars 

13055 H/S 6.54±0.57 



F 

N 
A 



6.62 
6.96 
6.81 



7.28 
6.01 
7.72 
6.27 
7.41 
6.25 



31205 H/S 8.25±0.92 
8.55 



F 

N 

A 
O 
T 

52271 H/S 
F 

N 

A 



7.09 
6.05 
6.61 

6.51 



7.66 
4.77 
7.93 
5.51 

7.2 



3.4±0.71 
, n , 4.67 
3 - 76 3.03 
q o q 4.33 
2.55 
4.26 



3.43 



2.77 



56731 H/S 7.07±0.68 
, 7.67 



F 

N 

A 
O 



6.89 
7.27 
7.09 



6.18 
8.16 
6.48 
7.80 
6.44 



18.78±0.54 
10.76±2.2 

8.95±2.4 

10.00±2.0 

-15.95±0.89 
-16.99±1.2 

-18.35±1.4 

-17.36±1.1 

-17.6±0.58 

-46.75±0.49 
-47.21±0.6 

-46.70±0.6 

-46.77±0.5 

-47.33±3.04 
-48.53±3.9 

-44.94±3.7 

-47.69±3.2 



-78.4±0.62 
-76.39±2.9 

-76.21±3.4 

-76.29±2.6 

6.22±0.84 
6.62±1.0 

7.47±1.1 

6.77±0.9 

6.7±0.49 

2.08±0.52 
1.99±0.6 

2.33±0.7 

1.87±0.6 

-26.48±0.5 
-26.26±2.2 

-29.18±2.5 

-27.24±2.0 



9.90±2.0 
11.95±2.1 
10.83±1.8 

4.26±1.1 

6.04±1.2 

4.72±1.0 
3.94±0.9 
5.22±0.63 

2.04±0.8 
2.17±0.9 
2.03±0.8 

7.91±1.3 

8.97±1.6 

7.88±1.2 
7.92±1.0 



86.3±13.4 
89.2±11.8 
87.1±10.8 

74.1±12.1 

86.6±8.8 

79.5±9.5 
74.8=1=11.9 
85±3.7 

69.4±26.1 
106.7±25.8 
76.1±22.9 

81.3±30.6 

69.2±29.1 

77.6±26.4 
73.7±24.3 



120=1=11 
117.5±14.4 

106.6±12.8 

113.1±11.7 

73.73±32 
98.5±20.4 

120.7±14.7 

109.5±16.2 

118±6.7 

57.5±25.9 
49.0±27.6 
49.3±22.4 

297.3±4.1 
298.1±19.8 

332.2±21.2 

308.8±17.6 
301.0 



281.1±15.0 
278.5±15.5 
280.5±12.6 

299.7±13.2 

295.9±8.9 

295.2±9.9 
313.8=1=12.7 
297±4.9 

30.6±24.9 
38.2±24.3 
35.2±20.3 

358.5=1=18.8 

334.3=1=18.1 

347.8±15.9 
351.3±12.6 



1.2±0.42 

1 50 L98 
LM 1.09 

2.24 



1.72 
1.59 



1.28 
2.01 
1.22 



0.44±0.15 

0.60°;^ 

!-00 IE 
71 1M 

u./i 47 

0.48±0.11 
0.98 



0.78 



0.61 



0.69±0.25 

54 °- 94 
U - M 0.26 

0.65 q 29 
1.01 



0.59 



0.30 



1.15 , 



1.1±0.37 
1.49 
0.86 
, ~, 1.62 
LIi 0.91 

1 11 141 

0.86 

1.1 ±0.2 



5.83±0.13 78 

8 09 10 - 71 
* w 5.62 



9.31 
8.61 



12.15 
6.75 
10.91 
6.44 



9.03±0.07 
13.74 f lf 

23.70 

2478 



16.69 



10.41 



7.64±0.09 

Am 11-19 
6.03 2 14 

7.39 

11.96 



6.79 



2.81 



6.07±0.08 
9.i 



7.36 
7.47 
7.04 



4.31 
10.53 
3.79 
9.13 
4.34 



38 


130.8 


46. 


.8 





40 


110.1 


37. 


.6 





80 


192.9 


85. 


.3 





79 










39 


61.8 


26. 


.8 





40 


58.6 


23. 


.9 





80 


96.3 


54. 


.7 






191 








101 








49 


42.8 


33.4 





52 


54.4 


47.0 





104 


108.1 


99.3 


4 


71 








34 


170.6 


25.9 





37 


146.3 


37.9 





74 


240.7 


77.1 






Table 2. (Continued.) 



HIP 



Sol 



VJ 

(mas) 



(mas/y) 



Us 
(mas/y) 



a 
(mas) 



i 

(°) 



(°) 



(°) 



ao/zu 
(AU) 



(km/s) 



N 



X X a 
HIP Orb (%) 



Strong Ba stars (cont.) 

103546* H/S 
F 

N 

A 

H/S 
F 

N 

A 



110108 



3.16±0.75 
4.20 
2.56 
4.17 
2.17 
4.13 
2.55 



3.28 
3.01 
3.24 



4.59±1.18 
7.55 



6.05 
5.65 



5.34 



4.84 
7.50 
4.26 
6.77 
4.21 



-2.89±0.82 
-3.08±0.9 

-2.63±1.1 

-3.01±0.8 

6.02±2.01 
13.29±6.9 

14.88±4.6 

10.16±4.0 



-20.06±0.67 
-19.77±0.7 

-20.76±0.9 

-19.97±0.7 

21.09±0.74 
18.75±1.9 

18.71±1.4 

19.77±1.2 



1.44±0.9 59.6±49.5 

1.99±1.0 80.4±35.1 

1.42±0.8 45.8±62.7 

7.16±4.8 92.4±16.2 

7.50±2.9 79.5±15.3 

5.30±2.6 85.1±19.3 



47.9±2.6 
141.4±52.9 

103.4±31.8 

132.9±82.6 

73.65±14 
151.1±29.4 

181.8±26.7 

172.3±28.4 



171.8±55.4 
234.6±35.8 
185.1±87.6 

118.3±15.7 
125.0±18.5 
128.2±21.3 



0.51±0.18 
0.91 



0.44 
0.66 
0.44 



0.13 
1.37 
0.24 
0.87 
0.15 



0.74±0.26 
2.48 
0.31 
2.45 
0.61 
1.87 
' 0.40 



1.18 j 
1.33 ; 
0.99 , 



7.77±0.11 
17.72 
0.45 
26.60 



7.37 , 



12.62 
6.08 



3.29 
16.93 
0.00 



5.1±0.06 
26.50 
3.12 
26.25 
5.84 
19.96 
3.95 



12.65 
13.96 
10.58 



70 



Remark: HIP 103546: Orbit not well defined 



35 


32.4 


28.8 


12 


35 


27.1 


21.8 


1 


72 


54.1 


49.8 


26 


25 








12 


13.3 


4.6 





13 


19.0 


8.4 


1 


26 


27.1 


18.9 
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Table 2. (Continued.) 



HIP 



Sol 



(mas) 



Ma 
(mas/y) 



H6 
(mas/y) 



a 
(mas) 



i 

o 



o 



0. 

O 



ao/zu 
(AU) 



(km/s) 



N 



X 

HIP 



X 

Orb 



(%) 



Tc-poor S stars 

17296 H/S 
F 



25092 



99312* 



6.07 



N 7.02 



A 
O 
T 

H/S 
F 

N 

A 

T 

H/S 
F 

N 

A 



6.27±0.63 
6.87 



6.29 



7.9 



5.37 
7.85 
6.27 
6.98 
5.66 



7.6 



1.32±0.99 
3.39 
0.88 
3.84 
1.14 
3.40 
1.13 
4.3 



1.73 ; 

2.09 
1.96' 



3.4 



2.7 



2.25±0.77 
3.78 8;2| 
9.29 



4.54 
4.18 



2.21 
7.98 
2.18 



-16.97±0.52 
-16.11±0.6 

-18.14±0.6 

-16.98±0.5 

-17.9±0.21 

48.18±0.79 
48.49±1.0 

48.85±1.0 

48.58±0.9 

49.8±0.84 

7.43±0.54 
7.44±0.6 

6.88±0.7 

7.28±0.6 



19.34±0.61 
19.24±0.7 

19.27±0.7 

19.35±0.6 

19.3±0.22 

-5.5±0.63 
-5.84±0.8 

-5.85±0.9 

-5.82±0.8 

-6.5±0.69 

1.1 1±0.52 
1.06±0.6 

1.23±0.7 

1.03±0.5 



1.72±0.7 

3.77±0.9 

2.30±0.6 
1.80±0.6 
4.50±0.29 

2.26±1.0 
1.47±1.1 
1.84±0.9 
2±1.1 

2.02±1.9 
2.96±2.1 
2.42±1.8 



85.6±25.3 

117.1±15.9 

105.6±18.7 
115.6±22.7 
114.±4.3 

97.5±22.4 
106.7±41.1 
100.2±26.2 
115±21 

113.9±41.8 
132.4±31.5 
122.5±32.0 



322±15 
335.7±26.4 

346.5±17.4 

334.3±19.0 
322 
361.±4.5 

232.51±33 
355.7±36.3 

322.0±61.9 

351.4±40.7 

306±24 

128.63±1.1 
80.3±60.0 

80.6±59.4 

82.4±49.6 



122.5±29.4 

194.0±16.8 

162.1 ±20.1 
166.8±24.5 
217±4.9 

302.4±25.9 
299.9±48.1 
300.6±29.8 
297±20 

253.4±42.4 
252.0±42.1 
255.1±34.3 



0.51±0.17 
, 0.44 
0.15 
0.74 
0.37 
0.52 
0.24 



0.28 
0.54 
0.37 



0.28±0.1 



0.57 



0.63 
0.51 



1.31 
0.71 
0.94 
0.7 



0.38 

2.27 

0.10 

2.46 

0.27 

1.3 

0.3 



8.47±0.2 
8.13 
2.42 
13.28 
4.97 
9.46 
3.63 



5.16 
8.77 
6.46 



0.53±0.18 
3.69 



7.888±0.276 
62.81 



22.07 
11.50 
15.75 



5.54 

38.68 

0.87 

41.81 

3.67 



0.37±0.13 22.3±0.09 



0.53 
0.65 
0.58 



2.28 
0.01 
2.30 
0.09 
1.94 
0.08 



14.57 
14.39 
14.56 



68.01 

0.17 

67.18 

0.76 

57.92 

0.97 



74 








36 


32.9 


22.7 





38 


71.7 


35.2 





78 


103.3 


83.9 


1 


154 








52 








26 


19.9 


16.4 


7 


26 


18.8 


17.4 


48 


52 


36.8 


34.1 


51 


72 








82 








40 


27.2 


25.8 


44 


42 


32.7 


29.9 


16 


84 


60.9 


59.7 


82 



Remark: 

HIP 99312: orbit not well defined 
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D. Pourbaix & A. Jorissen: Hipparcos parameters of Ba, CH and S stars 



HIP50805 with data from FAST 



HIP31205 with data from NDAC 



HIP1 16233 with data from FAST 




-ls.smas ^*=135.34459518deg. L8 - 8 ™» 

HIP 103546 with data from both consortia 




-9.2 mas ^*=224.77492779 deg. 

HIP105881 with data from NDAC 




297.37139926 deg. 




-7.3 ma., £*=96.33640079 deg. 

HIP62409 with data from NDAC 




-15.9 mas ^*=1 83.63073227 deg. «•»"■ 

HIP53763 with data from both consortia 



\ 








~Y / 





^*=123.38385363 deg. 




-123 mas ^*=345.54249700 deg. 

HIP1 16233 with data from NDAC 




-14.2 m« ^*=345.54249853 deg. 142 m » s 

HIP1 16233 with data from both consortia 




5*=345.54249767 deg. 



Fig. 1. The orbital arc on the plane tangent to the line of sight for some representative orbits among those listed in Table |[ 
The segments connect the computed position on the orbit to the great circle (not represented, perpendicular to the segment) 
corresponding to the observed position (Hipparcos measurements are one-dimensional). HIP 50805 is the only system for which 
the DMSA/O provides an orbital solution from scratch; HIP 53763 has a small parallax (about 2 mas), not very well determined 
since the orbital period is close to 1 yr; only the NDAC solution is acceptable for HIP 62409; the orbit of HIP 103546 is at limit of 
what can be extracted from the IAD; HIP 105881 is an example of an incomplete, albeit well determined, orbital arc; HIP 1 16233 
is one case where the NDAC and FAST solutions are rather different 



way to resolve that discrepancy is to assume that the orbital 
inclination is largely in error. 

The semi-major axis R — ao /w as derived from its astro- 
physical estimate (Eq. |l3[) is compared to its astrometric value 
in Fig. H and the two values are often consistent with each 
other. Although the value of an is likely to be affected by a pos- 



itive bias (i.e., a positive ao value is derived even when the data 
consist of pure noise, as clearly apparent from the astromet- 
ric ao/vo values listed in Table [j]), this bias does not markedly 
affects the retained solutions displayed on Fig. [| except for so- 
lutions with m < 3 mas, which all have R/R > 1. Solutions 
for larger zu values are almost equally distributed around unity. 



D. Pourbaix & A. Jorissen: Hipparcos parameters of Ba, CH and S stars 
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Fig. 2. Comparison of the astrometric (ordinate) and spectro- 
scopic (abscissa) determinations of the argument of periastron 



Fig. 3. Comparison of the semi-amplitude of the radial veloc- 
ity curve derived from the astrometry (using Eq. |l2|) and the 
spectroscopic value. The three panels show the results derived 
from FAST, NDAC and FAST+NDAC (from top to bottom re- 
spectively) 



3.3. Comparison TD vs IAD 

Among the 81 systems studied in this paper, 22 have TD and 
only 5 belong to the list of accepted orbits (HIP 17296, 25092, 
3 1205, 5080, and 105881). They hence can be used to check the 
consistency of the two orbit-determination methods. The astro- 
metric orbits derived from the TD are given in Table || as lines 
labeled with 'T' in column 2. The agreement between the re- 
sults obtained with the two data types is quite good (within la 
error bars). The number of TD measurements for those five sys- 
tems is up to twice as large the number of IAD measurements. 
As a consequence, the confidence intervals of the parameters 
is systematically narrower. One should however keep in mind 
that more precise does not imply more accurate. For instance, 
in the case of the mild Ba star HIP 105881, the 86 TD measure- 
ments yield ao/zu = 1.59 (with la limits of 1.45 and 1.74) as 
compared to 1.2 ± 0.35 estimated from Kepler's third law. The 
54 IAD measurements from FAST and NDAC combined yield 
a more accurate result of 0.99 (between 0.44 and 1.69). 

For the sake of completeness, we checked the astrophysical 
consistency of the orbit derived for the 22 systems for which 
TD are available. All systems but the five already accepted were 
rejected. In essence, this confirms that the astrometric content 
of the TD is basically equivalent to the IAD and that nothing 
new can come out from TD if IAD do not yield a reliable solu- 
tion. 



4. Comparison of parallaxes and proper motions derived 
from single- or double-star solutions 



It has been stressed (Wielen 1997; Wielen et al. 1999) that the 



proper motions provided by the Hipparcos catalogue may be 
systematically in error for long-period binaries (i.e. with P > 
3 yr), if those were not recognized as such by the reduction 
consortia. The orbital motion may then add to the actual proper 
motion, changing both its direction and modulus. The present 
sample offers a good opportunity to evaluate the impact of this 
effect. 

Figs. ^ and |§] compare the position angle and modulus, re- 
spectively, of the Hipparcos proper motion with those derived 
when account is made of the orbital motion as in the present 
work. It is clearly apparent that the direction of the proper mo- 
tion listed in the Hipparcos catalogue is correct, as expected, 
for orbits with periods less than 3 yr, corresponding to the du- 
ration of the Hipparcos mission. The Hipparcos values are the 
less accurate in the orbital-period range 3 to about 5 yr. The 
position angle quoted by the Hipparcos catalogue may then 
be off by several dozens degrees (a good illustration of this 
situation is offered by the the strong Ba star HIP 110108 in 
Table ||), though it differs generally by less than 3cr from the 
correct value. For longer orbital periods, the orbital motion be- 
comes negligible over the duration of the Hipparcos mission, 
and the proper motions are again rather well determined in the 
Hipparcos catalogue. The situation is almost identical for the 
the proper-motion modulus, except that the error bars now be- 
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Fig. 4. Uncertainty on the inclination j as a function of i. All 
points lying inside the region delineated by the oblique dashed 
lines have orbital inclinations consistent with i = 90°. Filled 
triangles correspond to orbital solutions derived from NDAC, 
filled squares to FAST and crosses to NDAC+FAST 



come very large for orbital periods longer than 3 yr. This sit- 
uation translates the fact that the proper motion and the semi- 
major axis are strongly correlated (resulting in a large formal 
uncertainty on the proper motion), because the two motions 
are difficult to disentangle when the Hipparcos data sampled 
only a fraction of the orbit. This situation is encountered e.g. 
for HIP 104785 and 104732 (whose astrometric orbits have 
not been retained precisely because of the large uncertainty on 
do introduced by the strong correlation with /i), and results in 
sometime large differences between the ao and \i values de- 
rived from the NDAC, FAST and NDAC+FAST data sets, since 
they are now very sensitive to the measurement errors. 

Fig. ^ compares the parallax listed in the Hipparcos cata- 
logue with that derived in the present work, as a function of 
the orbital period. It turns out that the agreement is good, ex- 
cept for systems with orbital periods close to 1 yr. For those 
systems, the parallax cannot be accurately derived since the or- 
bital motion and the parallactic motion are strongly entangled. 
This situation is encountered for HIP 17402, 53763 and 62827. 
Large error bars at periods different from 1 yr correspond to 
systems with parallaxes smaller than 3 mas. 

A few stars (HIP 8876, 29740, 29099, 32831, 36613) in 
our sample had negative parallaxes in the Hipparcos catalogue. 
The parallaxes we derive for these systems are slightly positive 
(w' was adopted to guarantee that property), but the associated 
error bar encompasses zero, so that these new parallaxes are 
just useless. 



Fig. 5. Correlation between the uncertainties on the inclina- 
tion i and the semi-major axis ao. The symbols have the same 
meaning as in Fig. f| 



5. Masses 



The test presented in Sect. [2.1.5| (Eq. |13|) to evaluate the astro- 
physical relevance of the orbital solution on the basis of ao /w 
is in essence based on an a priori knowledge of the masses, 
since Eq. O may in fact be rewritten as 



Ml 



(Mi + M 2 



p2 



The orbital solutions that appear to be statistically signif- 
icant on the basis of the F-test in almost all cases turned out 
to have a semi-major axis a^/w consistent with its expected 
value based on the estimated masses, the dwarf barium star 
HIP 60299 being however a notable exception. This agreement 
is well illustrated on Fig. |[ 

Therefore, the present analysis does not add much to our 
previous knowledge of the masses, especially since the astro- 
metric orbit only allows to eliminate i from the mass function 
but does not give access to the individual masses. 

It was originally hoped that the present astrometric results 
might be used to test the hypothesis which played a central 
role in the statistical analysis of the mass functions of CPRS 
( Jorissen et al. 1 998 ), namely that their Q distributions are quite 
peaked, because they host a white dwarf companion. However, 
the error bar on a^/ru (and hence on Q) is too large, even for 
the best determined barium-dwarf orbits (Table ^), to be able to 
draw a meaningful Q-distribution from the astrometric orbits. 
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Fig. 6. Comparison of the semi-major axis R = ao/w de- 
rived from the astrometry and from Eq. |l3| The upper panel 
displays the results from the FAST data, the middle panel the 
results from the NDAC data, and the lower panel the results 
from FAST and NDAC combined. The error bars assume that 
Mi, M2 and P are error- free 



6. Conclusions 

The availability of the Hipparcos Intermediate Astrometric 
Data and, when they exist, the Transit Data, allow to improve 
upon the existing Hipparcos catalogue. Once a new or revised 
spectroscopic binary orbit becomes available, it can be used to 
update the Hipparcos astrometric parameters and, sometime, to 
derive the corresponding astrometric orbit. In the latter case, 
the orbital inclination is obtained, thus allowing to improve our 
knowledge of the binary system. 

The optimization method used in the present paper yields 
a solution in all cases, and it is therefore important to perform 
consistency checks. One test which may be applied in all cases 
involves the comparison of the arguments of the periastron and 
of the semi-amplitude of the radial-velocity curve, as derived 
from either the spectroscopic orbit or the astrometric one. The 
comparison of the semi-amplitude of the radial-velocity curve 
involves however the inclination of the orbit which is not al- 
ways very well determined. 

In this context, chemically-peculiar red stars (CPRS) like 
barium, Tc-poor S and CH stars are interesting because they 
provide their own consistency checks. Indeed, beside the orbit 
which can be determined as for any other spectroscopic binary, 
the mass of both components may be guessed with some confi- 
dence. These mass estimates based on astrophysical considera- 
tions may then be compared to the mass of the system derived 
from the astrometric orbit using Kepler's third law. 



Fig. 7. Comparison of the position angle (expressed in radians) 
of the Hipparcos proper motion with that derived when account 
is made of the orbital motion (present work), as a function of 
the orbital period. The upper panel displays the results from the 
FAST data, the middle panel the results from the NDAC data, 
and the lower panel the results from FAST and NDAC com- 
bined. Error bars do not include the uncertainty on the Hippar- 
cos values 



From a sample of 8 1 CPRS spectroscopic systems whose 
orbits became available after completion of the Hipparcos cat- 
alogue, we have derived 23 reliable astrometric orbits. This 
shows that the number of Hipparcos entries for which an or- 
bital solution may be obtained is much larger than suggested 
by the number of existing entries in the DMSA/O. 

Updated astrometric parameters from this particular sam- 
ple have shown that the Hipparcos-catalogue parallaxes are not 
reliable for systems with orbital periods close to 1 yr. Similarly, 
the Hipparcos proper motions are not reliable for unrecognized 
binaries with periods in the range 3 to about 8 yr, as already 
suspected by |Wielen et al. (1999) . 

The orbit-determination methods based on the IAD and TD 
being completely different, systems for which the two sets of 
data are available can be used to assess the consistency of the 
solutions derived from both sets with different methods. It turns 
out that the two sets of results are generally in good agreement. 
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